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ABSTRACT 

Context. The solar mid-infrared metallic emission lines have already been observed and analyzed well, and the formation scenario of 
the Mgl 12 /jm lines has been known for more than a decade. Detections of stellar emission at 12 yum have, however, been limited to 
Mg I in very few objects. Previous modeling attempts have been made only for Procyon and two cool evolved stars, with unsatisfactory 
results for the latter. This prevents the lines' long predicted usage as probes of stellar magnetic fields. 

Aims. We want to explain our observed Mgl emission lines at 12 ^m in the K giants Pollux, Arcturus, and Aldebaran and at 18 jim 
in Pollux and Arcturus. We discuss our modeling of these lines and particularly how various aspects of the model atom affect the 
emergent line profiles. 

Methods. High-resolution observational spectra were obtained using TEXES at Gemini North and the IRTF. To produce synthetic 
line spectra, we employed standard one-dimensional, plane-parallel, non-LTE modeling for trace elements in cool stellar atmospheres. 
We computed model atmospheres with the MARCS code, applied a comprehensive magnesium model atom, and used the radiative 
transfer code MULTI to solve for the magnesium occupation numbers in statistical equilibrium. 

Results. The Mg I emission lines at 12 /urn in the K giants are stronger than in the dwarfs observed so far. We present the first observed 
stellar emission lines from Mg I at 1 8 jim and from Al I, Si I, and presumably Ca I at 12 /urn. We successfully reproduce the observed 
Mg I emission lines simultaneously in the giants and in the Sun, but show how the computed line profiles depend critically on atomic 
data input and how the inclusion of energy levels with n > 10 and collisions with neutral hydrogen are necessary to obtain reasonable 
fits. 

Key words. Infrared: stars, Stars: atmospheres, Stars: individual (Pollux, Arcturus, Aldebaran), Stars: late-type, Line: formation 



1. Introduction 

Metallic solar emission li n es aro und 12/zm were first iden- 
tified by Chang & Noyesl (1 19831). with the most prominent 
lines originating from transition^ 7/ — > 6h (12.32yum) and 
Ih — > 6g (12.22 jjm) between Rydberg states of neutral mag- 
nesium. Additional Rydberg emission li nes from All, Si I, and 
tentatively C al were identified as well (IChang&Noveslll983t 
IChangl [19841). The Mgl line form ation scenario r emain ed un- 
clear until IChanget all d!991l) and lCarlsson et all d 19921) . here- 
after C92, in two independent studies reproduced the emission 
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features by employing standard plane-parallel numerical radia- 
tive transfer with a detailed atomic model and a reliable solar 
atmosphere. They confirmed an origin below the atmospheric 
temperature minimum, refuted a chromospheric line contribu- 
tion, and established a non-LTE (LTE: Local Thermodynamical 
Equilibrium) formation scenario. The solar lines have subse- 
quently been used in, e.g., Mgl statistical equilibrium analy- 
sis by Zhao et al. (1998). C92 proposed a general Rydberg line 
formation mechanism for the highly excited metal lines, which 
implied that all visible metal emission lines in the solar spec- 
trum around 12/zm originated in the photosphere. A detailed 
non-LTE modeling of t he Al I emission has been carried out by 
Bau mueller & Gehrenl (119961) . where they confirmed this mech- 
anism. 

The lack of suitable spectrometers and the low stellar flux in the 
mid-infrared have in the past made high-resolution spectroscopy 
in this wavelength re gion possib l e only for the Sun and a few lu- 
minous nearby stars. iRvde et al.l(l2004l) observed the 12 /mi Mg I 
emission features in Procyon, and successfully reproduced the 
line profiles by employing th e same modeling technique as C92. 
lUitenbroek & Novesl (I 19961) observed and modeled the evolved 
stars Arcturus (a Boo) and Betelguese (a Ori). Using the same 
model atom as C92, they were unable to fit the line profiles 
of the 7/ — > 6h Mg I transition, which appeared both in emis- 
sion (Arcturus) and absorption (Betelguese). Their observational 
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sample also included five M giants and supergiants , in wh ich the 
line appeared in absorption. However, Rvd e et all (12006) inves- 
tigated water vapor lines for Betelguese in the same spectral re- 
gion, and found a water line that coincided with the wavelength 
of the Mgl 12.32/im line. The group successfully modeled the 
water line, without considering the Mg I blend (which we pre- 
dict to be very weak, see Sect. |6.31 >. This may explain the sample 
of observed M star absorption at 12.32 yum, since water vapor is 
expected in these stars, whereas the Mg I emission line contribu- 
tion should be minor. 

A well known potential use for the Mg I lines is as probes of 
magnetic fields, which play a fundamental role in the underly- 
ing physics of a cool stellar atmosphere. Zeeman line-splitting 
from an external magnetic field increases quadratically with 
wavelength, while the Doppler broadening only has a linear de- 
pendence. Thus a line's sensitivity to magnetic fields becomes 
higher at longer wavelengths. The splitting of the solar emis- 
sion lines was pointed out early and has been extensively ana- 
lyzed. We have performed observations of the magnetically ac- 
tive dwarf e Eridani. These will be reported on in a forthcoming 
paper (Richter et al., in preparation), hence we defer further dis- 
cussions about stellar disk-averaged magnetic fields until then. 
Prior to (stellar) diagnostic applications, however, we should 
make sure that we are able to model and understand these lines 
in a range of stars. So far, as mentioned above, modeling at- 
tempts for evolved stars have been unsuccessful. We address this 
issue here by analyzing high-resolution observational spectra, 
which show strong Mg I emission lines in the three giants Pollux 
(KOIII), Arcturus (K1.5III), and Aldebaran (K5III). We model 
and analyze simultaneously the three K giants and the Sun, with 
particular emphasis on influences from atomic data, and discuss 
why previous modeling attempts have not succeeded. The orga- 
nization of the paper is as follows; in Sect. [2] we describe the 
observations. In Sect. [3] we review some concepts about the for- 
mation of the infrared Mg I emission lines and in Sect. [4] we 
describe our modeling procedure. Results are presented in Sect. 
[5] and we discuss them and give our conclusions in Sect. [6] 

2. Observations 

The observations were ma de with TEXES, th e Texas Echelon- 
cross-echelle Spectrograph. lLacv et al.ld2002l) . TEXES provides 
high spectral resolution in the mid-infrared and is available as a 
visiting instrument at both Gemini North and at IRTF, the Infra- 
Red Telescope Facility. The Pollux (B Gem) observations come 
from the November 2006 observing campaign at Gemini North. 
The Arcturus and Aldebaran {a Tau) observations were done 
over many years at the IRTF. In most cases, the observations 
were primarily intended for flux calibration or focus tests and 
not to study the stars themselves. 

When observing stars with TEXES, we nod the source along 
the slit, typically every 10 seconds, to remove sky and telescope 
background. Before each set of 8 to 16 nod pairs, we observe a 
calibration sequence that includes an ambient temperature black- 
body and an observation of blank sky emission. The difference 
of blackbody minus sky serves as a first order telluric correc- 
tion and flatfield. Where possible, a featureless continuum object 
with emission stronger than the target is also observed to further 
correct for telluric features and flatfielding. The largest asteroids 
work very well for this purpose, as does Sirius (a CMa) with re- 
spect to Pollux. 

At the frequencies of the mid-infrared Mg I emission lines, the 
spectral orders from the TEXES high-resolution echelon grat- 
ing are larger than the 256 2 pixel detector array. This results in 



slight gaps in the spectral coverage. For the Pollux observations, 
which were the final observations before sunrise, we observed in 
two settings and adjusted the tilt of the collimator mirror feeding 
the echelon grating for the second setting. This shifts the spectral 
orders in the dispersion direction. By combining the data from 
these separate observational settings, we were able to fill in the 
gaps in the spectral orders. The Arcturus and Aldebaran data 
were constructed from many separate observing settings and no 
particular efforts were made to fill in the gaps. 
Data reduction w as done using a custom FORTRAN pipeline 
(lLacv et al.ll2002l) . The pipeline corrects for spikes and optical 
distortions in the instrument, allows the user to set the wave- 
length scale based on telluric atmospheric features, flatfields the 
data, differences nod pairs to remove the background emission, 
and then combines the resulting differences. Finally it extracts 
a spectrum based on the spatial information within the two- 
dimensional echellogram. The pipeline also provides a fairly ac- 
curate estimate of the relative noise in each pixel. 
To combine data from separate observations, we first established 
a common wavelength scale. We corrected each spectrum for the 
Earth's motion at the time of the observation and then interpo- 
lated the data onto the common scale. We used a fourth-order 
polynomial derived from line free regions to normalize each 
spectral order. We determined the signal-to-noise (S /N) for the 
normalized spectrum via a Gaussian fit to pixel values and used 
the relative noise estimate established during pipeline reduction 
to assign a weight for each spectral pixel. When combining data, 
we choose to weight by the signal-to-noise squared, which ef- 
fectively means weighting by successful observing time. 
Observations of low pressure gas cells near 13.7 yum at the 
November 2006 run indicate that the instrumental profile for the 
12 /mi observations of Pollux has a Gaussian core with a FWHM 
~ 3.0 km s , corresponding to a spectral resolution R ~ 10 . As 
the Arcturus and Aldebaran spectra combined data from four 
different runs and possible errors from the combinations may 
be significant, we were unable to make a reliable measurement 
of the instrumental profile in this region for these stars. At 18 
/mi, similar measurements indicate that these observations have 
a Gaussian instrumental profile with a FWHM ~ 4.5 km s . In 
Sect. [5] we display our observed data re-binned to approximately 
the spectral resolution, except for Fig. [5] where the pixel scale is 
used. Signal-to-noise ratios in the spectra vary but are generally 
high, reaching S /N ~ 450 per pixel for Pollux and ~ 300 for 
Arcturus and Aldebaran, in regions around the 12.22 /mi line. At 
18.83 fim, the ratio is S /N ~ 40. 

3. Departure coefficient ratios 

Before proceeding to a modeling description, we briefly review 
some important concepts about the formation of the 12 //m emis- 
sion lines. In the following we use the departure coefficients 
b/ = rij/n*, where n,- is the actual number density (not to be con- 
fused with the principal quantum number n) of energy level i and 
n* the corresponding LTE population, as calculated from the to- 
tal magnesium abundance using the complete Saha-Boltzmann 
relations. In a spectral line, a departure coefficient ratio which 
differs from unity, b\jb u + 1, at line-forming depths causes a de- 
viation of the line source function, S [, from the Planck function, 
B v , which affects the emergent intensity: 

cl „hv/kT_i 

— = (1) 

B v bilb u x e hv l kT - 1 

For a characteristic wavelength A - 12.3 /mi, and temperature 
T = 5000 K, we get e hy/kT ~ 1.26, and may directly from Eq.Q] 
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realize that already a small deviation from unity in the departure 
coefficient ratio causes a significant change in the line source 
function. The physical reason for this is the increasing impor- 
tance of stimulated emission in the infrared. If b„/bi > 1 and 
increases outwards in the atmosphere, we may get a rising total 
source function and a line profile appearing in emission despite 
an outwards decreasing temperature structure. Such departure 
coefficient divergence occurs between highly excited Rydberg 
levels in the outer layers of the modeled stellar photospheres 
considered in this study, and is the reason for the modeled emis- 
sion lines. 

Departure coefficient ratios that deviate from unity are set up by 
three-body recombination from the Mg II ground state and a 'de- 
excitation ladder' that preferably takes An, Al = — 1 downward 
steps (see Fig.[T]for an illustration). In the solar case, all Mgl 
Rydberg levels are strongly collisionally coupled to each other 
and to the Mgll ground state. The main effects that drive the 
line source function out of LTE come from lines elsewhere in the 
term diagram, primarily lines between levels with intermediate 
excitation energies, which are opticall y thin in the outer atmo - 
sphere and experience photon losses dRutten & Carlssoral 19941) . 
These levels impose a lower limit to the Rydberg state deexcita- 
tion ladder. The number densities of the Rydberg energy levels 
adjust to the upper and lower limits, and a radiative-collisional 
population flow occurs. It was shown, for the solar case by 
C92, how a high probability for An, A/ = — 1 downward tran- 
sitions is necessary for the Rydberg state ladder to be efficient 
and that these transitions dominate only if the collisional cou- 
pling is strong in the uppermost Mg I levels. It was also shown 
how this high probability arises from the regular character of 
the collisional cross-sections of transitions between highly ex- 
cited levels. We thus remind of the remark in C92 that for a cor- 
rect description of the ladder flow between highly excited lev- 
els, it is more important to have a consistent set of collisional 
data, than to have the most accurate cross-sections for a few 
transitions. These considerations are important to keep in mind 
when we later discuss our extension of the model atom and the 
influence from collisions with neutral hydrogen. A more com- 
prehensive description of non-LTE effects throughout the Mg I 
term diagram that affect the solar 12 /urn lines can be found in 
iRutten & Carlssonl d!994l) . 

4. Modeling 

To produce synthetic line spectra, we employ standard one- 
dimensional, plane-parallel, non-LTE modeling for trace ele- 
ments in cool stellar atmosphere s. We generate mode l atmo- 
spheres from the MARCS code dGustafsson etalj|1975l 120081) . 
adopt a comprehensive magn esium model atom, an d use the ra- 
diative transfer code MULTI dCarlssonll 19861 1 19921) to solve for 
the magnesium occupation numbers in statistical equilibrium, 
while holding the structure of the atmosphere fixed. 

4.1. Model atmospheres and stellar parameters 

The MARCS hydrostatic, plane-parallel models are computed 
on the assumptions of LTE, chemical equilibrium, homogene- 
ity, and the conservation of the total flux (radiative plus convec- 
tive; the latter treated using the mixing-length theory). No chro- 
mospheric temperature rise is invoked but, as shown by C92, 
omitting a chromosphere has a negligible impact on the solar 
Mgl 12 jum trans itions. The findings in this work and that by 
Rvde et all d2004l) suggest that this holds true also for other in- 
vestigated stars. Apart from the temperature and density stratifi- 
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Fig. 1. Illustration of the 'Rydberg ladder' (see text), using a 
selected part of the Mg I term diagram. Five levels about the 
Ih — > 6g transition are labeled with corresponding quantum 
numbers nl. Dots mark energy levels with n > 9 and I > n — 3. 
Solid lines between levels show transitions with An, Al = —I, 
dashed show two alternative transitions (An = 2, Al = -1). 
Dashed-dotted illustrate recombination from the Mg II contin- 
uum. 



cations, a detailed MARCS radiation field was generated using 
opacity samplings including millions of lines. A sampled ver- 
sion of this radiation field was used when MULTI calculated the 
photoionization rates, in order to properly account for the line- 
blocking effect. We discuss some variations to our procedure in 
Sect.|6] 

The stellar parameters we used for Pollux were T e g = 
4 865 K, log g = 2.75 (cgs), a solar metallicity (as given by 
iGrevesse et alj 2007), and a depth independent 'microturbu- 
lence' f = 1.5 km s _1 , all based on a spectral an a lysis o f optical 
iron and calcium lines made bv lDrake& Smith! d!99 11) . For the 
parameters of Aldebaran we adopted T e g = 3 900K, logg = 
1.5 (cgs), a metallicit>0 [ M/H] = -0.25 and f = Ukms" 1 . 
These are from primarily iDecin et ail (120031) but considering 
also sources accessible at the SIMBAD astronomical database. 
Finally for Arcturus we used T e s = 4 280 K, log g = 1.5 (cgs), 
[M/H]= -0.50, and £ = 1.7kms~' A discuss i on of these 
Arcturus parameters can be found in iRvde et al.l d2002b . The 
stars are all nearby and well-studied objects, and their param- 
eters should be fairly accurate. Model grids show that the ef- 
fects on the lines from (reasonable) variations in log g or T e g are 
smaller than effects from, e.g., atomic input data, which will be 
investigated in the following sections. 

We convolved our computed intrinsic line profiles with the 
instrumental profile, the projected rotational velocity (v sin i), 
and the 'macroturbulence' (none of which affect the line 
stren gth but only the p r ofile s hape). We adopted v sin/ values 
from Smith & Dominv ( 1979), which for Pollux, Arcturus and 
Aldebaran are, respectively, v sin i = 0.8, 2.7, 2.7 km s _1 . As we 
were unable to obtain a fair estimate of the instrumental pro- 
file for the Arcturus and Aldebaran spectra around 12 /mi, we 
choose first to assign an isotropic Gaussian shape with charac- 



2 Here we do not specify individual metal abundances, but [M/H] 
is taken as the characteristic metallicity where as usual [A/B] = 
log(n A /n B ), - log(n A /n B ) . 
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teristic Doppler velocity v m for the combined effect of the in- 
strumental profile and the macroturbulence. For Pollux, where 
the instrumental profile could be separated out, we obtained 
Vmacro ~ 3.3 km s _1 . However, it became clear that the modeled 
line wings of the K giants better fitted the ob servations w hen as- 
suming a radial-tangential Gaussian shape dGravlll976l) for the 
macroturbulence. Therefore we decided to assign the Pollux in- 
strumental profile for all three stars (the exact values are not so 
significant since the macroturbulence is the dominating exter- 
nal line broadening), and adopted a radial-tangential macrotur- 
bulence v m ji-T = 5.5, 6.0, 5.5kms _I to fit the observed line- 
widths. Our values are ~2kms 1 higher tha n those measured 
from F ourier analysis in optical spectra by ISmith & Dominvl 
(119791) . 

4.2. The model atom 

Our Mgl model atom is essentially an enlarged and slightly 
modified version of the one compiled by C92, and a full de- 
scription can be found there. The original model atom has 
also been u sed in the analysis of solar magnetic fields (e.g., 
Bruls et all 1995), for the Mgl 12/mi flux profiles of Procyon 
(Rvde et aLll2004l), and in a previous attempt to model giants 
dUitenbroek & Novesll 19961) . In short, the atom is complete with 
all allowed transitions up to principal quantum number n — 9 and 
includes the ground state of Mg EL We now describe changes and 
tests we have made. 



4.2.1 . Enlargement of the model atom 

iLemke & Holwegerl dl987h . who also pointed out a possible pho- 
tospheric line-origin through a rising line source function, made 
a statistical equilibrium investigation for the Sun but were unable 
to reproduce the Mgl 12/im emission due to a combination of 
their adopted collisional data treatment and an inadequate model 
atom. Their exclusion of levels higher than n = 7 resulted in an 
incorrect description of the replenishment of the Rydberg levels 
from the ion state. C92 experimented with the n — 8, 9 levels and 
confirmed that these were necessary to have top-levels that were 
fully dominated by collisions and to obtain sufficient departure 
coefficient differences in line-forming layers to match the ob- 
served solar emission features. In this work we have extended 
the model atom to include levels with n > 10 to investigate if the 
departure coefficient differences are further enhanced. 
The model atom was first enlarged to include all energy levels 
and allowed transitions with n = 10. To ensure homogeneity 
throughout the model atom, all new atomic data were calculated 
using the same formalisms as those employed by C92. The only 
exception was absorption oscillator strengths for transitions with 
/ < 3 for w hich data was drawn from the opacity project (OP) 
TOP-BAS E dCunto & Mendozal[l992l) . since the tabulation used 
by C92 dMoccia & Spizzoll 19881) only extend to n = 9. The en- 
largement caused an upward shift in the Mg I departure coeffi- 
cients, and the effect became more pronounced as n increased; 
thus producing larger departure coefficient differences between 
adjacent levels. The same effect was seen in the Mgl statisti- 
cal equilibrium for all our template atmospheres. To investigate 
the influence of the low I levels, a test-run was also made where 
only n — 10 levels and transitions with I > 4 were included. This 
model atom and the complete n — 10 atom produced almost in- 
distinguishable results. 

After this initial enlargement, atomic models were constructed 
step-wise, including higher principal quantum numbers. The en- 
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Fig. 2. Solar departure coefficients for Rydberg state energy lev- 
els with n > 5, as functions of the continuum optical depth at 
500 nm. b$ is indicated in the figure, and an increment of one 
follows upwards in the figure. The displayed coefficients com- 
bine all b„ i coefficients with the same n and I > 4. The Mg II 
ground state is labeled with dots. 



hancement continued until the atom's uppermost n level and 
the second uppermost were Boltzmann populated with respect 
to each other at all atmospheric layers (i.e., b top ^\jb top = 1). 
This criterion was met for all atmospheric models when reach- 
ing n = 15, illustrated for the solar case in Fig. [2] (where our 
model actually meets the criterion already at n — 12). Sensitivity 
tests verified that no differences in results occurred when adding 
the final n — 15 top-levels. The main difference between the 
solar departure coefficients (Fig. p2]l and those of the K giants 
is that the latter have Mgll ground states that are more over- 
populated relative to LTE in their outer atmospheres. This is 
mainly because Mg I and Mg II are competing ionization states 
in these cooler atmospheres, so that the ion ground state becomes 
more sensitive to deviations from LTE in Mg I population den- 
sities. In Arcturus and Aldebaran, the overionization is further 
enhanced by the lower metallicity, which reduces the important 
line-blocking effect. 

As the influence from levels with low I was negligible already 
for n = 10, we have confined the enlargement to levels and tran- 
sitions with / > 4. 

4.2.2. Collisional data 

Rvde et al.l d2004l) compared radiative bound-bound and bound- 
free data from our model atom with OP data and found an overall 
good agreement. We have therefore restricted our analysis here 
to some aspects of the collisional data, which usually introduce 
the largest uncertainty in the model atom, due to a large number 
of poorly known cross-sections. 

4.2.3. Collisions with neutral hydrogen 

The role of collisions with neutral hydrogen in cool stellar atmo- 
spheres has long been a subject for debate. Despite small cross- 
sections (as compared to electron impacts), one may expect them 
to contribute significantly to collisional rates due to large rtHi/n e 
ratios. In the outer parts of the model atmospheres in this study, 
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this ratio ranges from about 10 4 in the Sun to 10 5 in the cooler 
and more metal-poor Arcturus. No collisions with neutral hydro- 
gen were considered in the original model atom. 
When inelastic collisions with neutral hydrogen are explic- 
itly included in non-LTE calc ulations , a st andard proce- 
dure is to adopt the recipe of iDrawi 3 (Qj69), as given by 
ISteen bock & Holwegerl d 19841). The Dr awin formula has often 
been criticized. ISteenbock & Holwegerl state an accuracy of an 
order of magnitude, but a rather common remark is that the 
recipe may overestimate the cross-secti ons with as much as one 
to six orders of magnitudes (see, e.g., lAsplundl 2005, and ref- 
erences therein). Unfortunately, more reliable cross-sections are 
scarce, especially for non-LTE calculations that require data for 
a large set of transitions. A customary way around this shortage 
is to adopt a scaling factor Sh to the Drawin formula, calibrated 
on solar or stellar observations. 

In studies concernin g the solar Mg 112 ^m lines, the Drawin for- 
mula w as adopted bv lLemke & Holwegerl d 1 987b and lZhao et alj 
d!998h . The latter group scale their values with a factor 
that decreases exponentially with increasing excitation energy. 
Consequently, they apply Sh = 3 x 10~ 10 for the 12 fim transi- 
tions, which give them essentially the same result as if neglecting 
hydrogen collisions. Recall also that the former group was un- 
successful in producing an emission line core in their statistical 
equilibrium analysis. 

Here we have estimated the collisional rates due to neutral hy- 
drogen impacts using the Drawin formula. When introduced 
without scaling factor, our models reveal solar 12 //m lines 
in pure absorption and this case will therefore not be consid- 
ered. By calibrating on the solar observations, we have adopted 
Sh - 10~ 3 in all computations involving collisions with neu- 
tral hydrogen. However, in view of the existing uncertainties, 
we present results both from including these collisions for all 
radiatively allowed bound-bound (for which we have oscillator 
strengths) and all bound-free transitions, and excluding them. 

4.2.4. Collisional excitation from electrons 

The bulk of the collisional cross-sections from electron impacts 
for radiatively allowed transitions are calculated, as was done i n 
C92, using the imp act parameter approximation dSeatonlll962h . 
Mashonkina ( 1996) showed that, overall, this approach predicts 
significantly smalle r cross-sections than th e alternative semi- 
empirical formula of lyan Regemorteri (1 1 9621) (applied for th e so- 
lar Mg I lines in, e.g. JZhao et all 1998V lAvrett et all d 19941) used 
both formalisms when modeling the solar Mg I lines and con- 
cluded that using Ivan Regemorteri gave somewhat weaker emis- 
sion. In C92, the impact parameter approximation is claimed to 
give a consistent set of rates accurate to within a factor of two 
for transitions between closely spaced levels. 
Note that the above mentioned formalisms relate the collisional 
cross-section to the oscillator strength and may therefore not 
be applied to radiatively forbidden transitions. The 'forbidden' 
cross-sections are here set to a multiplying factor times that of 
the closest allowed ( see C92). The orig inal choice (C92) for this 
factor was 0.05, but iBruls et al.l (Q995) discovered some errors 
regarding a few oscillator strengths, accounted for here as well, 
and the factor was revised to 0.3 in order to reprodu ce the pre- 
vious results. We also adopt 0.3, which was used by lRyde et al.l 
(2004) as well . A similar t r eatme nt for solar Mgl analysis has 
been used by iMauas et al.l (1 1988]) who assumed 0. 1, the same 
value as estimated in lAllenl (Il973l) . ISigut & Lesterl (j.1996) also 
adopted 0. 1 for Mg II (in work where they, for B type stellar pho- 
tospheres, predicted the corresponding Rydberg emission lines 



for Mg II), which they found to be in rough agreement with a 
few more rigorously calculated rates from low excitation tran- 
sitions. We have tested using the enlarged model atom without 
collisions with neutral hydrogen and concluded that by a raise to 
0.7 times the cross-section of the nearest allowed transition we 
are able to reproduce the observed solar lines, but that the mod- 
eled emission in the K giants remain far lower than the observed. 
An alternative approach for the radiatively forbidden transi- 
tions, m£\iedjn,ej' Mgl (and II) non-LTE abundance anal- 
yses riZhao et al.lll998tlPrzybilla et alJl200lHGehren et al.ll2004l: 
iMashonkina et al. 2008), is to set a constant collisional strength 
Q = 1. Overall, this gives considerably lower collisional rates. 
The collisional strength Q for collisions with electrons is related 
to the Maxwellian averaged downward collisional rate C,, [s~'] 
via: 

C fi = 8.63 x lO- 6 £lT; l/2 gJ l n e (2) 

where T e is the electron temperature, gj the statistical weight of 
the upper level, and n e the number density of free electrons. The 
upward rate, Cy, then follows from the principle of detailed bal- 
ance. We have also tested assuming a constant Q = 1 on the 
solar and Arcturus model, and verified that the modeled emis- 
sion for Arcturus still is inadequate, remaining on the same low 
level as that displayed in Fig. [3] when collisions with neutral 
hydrogen are not included. The solar emission increases by this 
approach, but not enough to drastically change the suitable Sh 
factor. Overall, our experiments with the radiatively forbidden 
transitions tell us that the solar emission lines are somewhat sen- 
sitive to these rates, whereas the models of the K giants are much 
less responsive. 



4.2.5. l-changing collisions 

No explicit calculations of collisional transitions of type 
«,/—>«,/' with high orbital quantum numbers (/' > 4) were made 
in the original model atom, but rather it was assumed that these 
collisional rates were high enough to ensure strong coupling be- 
tween closely spaced levels. This was established by setting Q = 
1 5 if the difference in effective principal quantum number (n - 6 
where 6 is the quantum defect) in the transition was less than 0. 1 . 
The high value of Q. causes very high collisional rates, which es- 
sentially force all levels with the same n and / > 4 to share a com- 
mon departure coefficient (i.e., to be Boltzmann populated with 
respect to each other), in agreement with the proposed assump- 
tion. This, however, only holds for all relevant atmospheric lay- 
ers in the Sun, whereas in the more diluted model atmospheres 
of the K giants, we find deviations between departure coeffi- 
cients in outer layers (for Arcturus outside log T500 ~ -2) for 
levels with relatively low orbital quantum numbers, when set- 
ting CI = 10 5 . Apparently, such a value does not suffice for the 
K giants and we need to either raise the factor or explicitly esti- 
mate the rates. Note that C92 based the assumption of common 
departure coefficients mostly on the large cross-section s for l- 
chang ing collisions with neutral hydrogen calculated bv lOmonll 
(Il977l) . cross-sectio ns later shown to be over-estimated b y an or- 
der of magnitude (Hoan g-Binh & van R egemorter 19951). 
In this work, we have calculated explicit electron/ion collisional 
rates for transitions with Al = ±l,An = and 1,1' > 3 
by using a cut-off at lar ge impact parameters, as outlined by 
iPengellv & Seatonl dl964l) . We find that the radial cut-off for 
transitions with low / is set by the non-degeneracy of the en- 
ergy levels. For ion rates, this happens inside the radius where 
the strong interaction dominates and hence the impact parame- 
ter approximation is not expected to be reliable. However, for 
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Fig. 3. Observed and modeled line profiles in Arcturus and 
Aldebaran for the 12.22 fim line, plotted on a velocity scale. 
Labels as indicated in the figure; where C92 denotes the orig- 
inal model atom, n=15;noH the extended excluding collisions 
with neutral hydrogen and n=15;H the one including such. The 
filled dots denote the observed data. 



electron impacts the cut-off is at larger distances and since elec- 
tron rates dominate over ion rates for these transitions, our ap- 
proach should provide rates accurate to at least an order of 
magnitude also in the low I range. Ion rates surpass electron 
rates from I ~ 9 and higher, where the cut-off is well within the 
weak-interaction limit. We assume only singly ionized elements 
(n wn = n e ) with the largest electron/ion donor being magnesium, 
providing ~ 40 % of the total electron/ion pool in the relevant at- 
mospheric layers in our MARCS models. 
The calculated rate s are in good agreement with the el ectron 
rates tabulated in Hoang-Binh & van Regemorter (1994), who 
exclusively considered Z-changing collisions for the n — 6,7 lev- 
els. Our rates are higher than those inferred from Q = 10 s (by 
typically a factor of ~ 4 for, e.g., n = 7), hence no changes in 
results occur for the solar atmosphere. For the K giants, small 
differences between departure coefficients with the same n and 
/ > 4 still exist but effects from using explicit rates are small. 
When forcing common departure coefficients for all equal n lev- 
els with I > 4 (by drastically increasing our computed rates), we 
still obtain normalized emission peaks for the K giants that differ 
only by a few percentage points. 



5. Results 

5.1. Emission lines at 12 /dm 

We plot observed and computed line profiles for the Mg 112 fim 
lines in the Sun, Pollux, Arcturus, and Aldebaran in Figs. [3] |4] 
and [5] In Arcturus and Aldebaran, the 12.32 fim line is b lended 
with a water vapor absorption line (see lRyde et al.l ((2006) for an 
identification), which in the latter star is so influential that we 
choose to exclude the Aldebaran 12.32 fim line from the analy- 
sis. The observed emission lines from the K giants are stronger 
than the solar lines. In addition to the Mg I lines, we also identify 
emission lines from Si I, Al I, and Ca I in the observed spectra of 
Pollux and Arcturus (three of the lines are displayed in Fig. [5]), 
all identified as n — 7 — > 6 transitions with high orbital quantum 
numbers. The flux maximum in the observed and normalized 



Fig. 4. Observed and modeled line profiles in Pollux and for 
three positions on the solar disk (indicated in the figure) for the 
12.22 fim line, plotted on a velocity scale. Labels as in Fig. [3] 



spectra (and the FWHM for the Mg I li nes) are given in Tables 
Q]and[2] Mo re line data can be found in Chang & Novesl (l983) 
and IChangl d!984h . We have also added results from modeled 
Mg I flux profiles for the Sun (from a disk integration over solar 
intensity profiles in a model that reproduces the observations) 
in Table [T] to enable a fair comparison between solar and stel- 
lar observations. This illustrates that the K giants have stronger 
emission than the Sun. 

The line-center average depth of formation in the modeled Mg I 
12 fim lines is, for the Sun and Pollux, in atmospheric layers 
slightly below logTsoo ~ -3, with the weaker 12.22 fim line 
shifted approximately 0.2 dex toward the inner photosphere. In 
Arcturus and Aldebaran, the line formation takes place deeper 
inside the atmosphere. The average depth of formation for the 
line-center in the 12.22 fim line in Arcturus and Aldebaran is 
logTsoo = —1.8 and -1.6 respectively. This is partly because 
of the lower amount of H ff opacity in these atmospheres (due 
mainly to lower electron abundances), which shifts the contin- 
uum formation to about logTsoo ~ -0.8, as compared with 
logTsoo 1.2 in Pollux. 

The extension of the model atom has a significant impact on 
the synthetic line spectra, with computed intensity/flux profiles 
being much stronger when using the enlarged model atom. We 
note the failure of the smaller atom to reproduce the observed 
emission for the K giants, whereas it provides a good match for 
the solar lines, in agreement with previous studies. The larger 
model atom without collisions with neutral hydrogen predicts 
emission lines well below the observed level for Arcturus and 
Aldebaran, in contrast to the Sun where the modeled lines now 
are too strong. However, when including collisions with neutral 
hydrogen (as described in Sect. 14.2.3b the models reproduce the 
observed emission in all cases. These different responses to the 
'added' collisions demonstrate the complexity of the Rydberg 
state deexcitation ladder, and are further discussed in Sect. 16.21 
Our models predict narrow absorption troughs in the Pollux 
lines, only matched by observations in the red wing of the 12.22 
fim line. However, due to uncertainties in the observed normal- 
ized spectra imposed by, e.g., the continuum setting, we are not 
able to draw any firm conclusions from the absence of absorp- 
tion troughs. A discussion about shifts in the s olar absorptio n 
troughs, visible in Fig. [5] can be found in, e.g., Chang (1994). 
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Fig. 5. Observed spectra around the 
12.3 fim region, shifted to the solar 
frame, for Pollux, Arcturus, and for 
two positions on the solar disk (in- 
dicated in the figure). Solar observa - 
tions are from lBrault & Novesl dl983l> . 
Models for the Mgl 12.32 yum line are 
labeled with dots, and use the n=15;H 
model atom (see Fig.O. The solar ab- 
sorption features are telluric lines, and 
the missing parts in the Arcturus spec- 
trum are due to gaps between spectral 
orders. Visible emission lines as indi- 
cated in the figure. The Arcturus ab- 
sorption line at 12.32 yum is a water va- 
por line. 



We also note how the line-wings in Arcturus and Aldebaran are 
too broad to be fitted by an isotropic Gaussian, and require a 
radial-tangential macroturbulence (see Sect. 14. lb . 

5.2. Mg I emission lines at 18fim 

We also observed the 8/i — > lg Mgl transition at 18.83 yum in 
Pollux and Arcturus, and present here the first stellar observa- 
tions of this line. The emission is high here as well, see Fig. [6] 
and Table [TJ Our synthetic line spectra reproduce the observed 
emission also for this line, which suggest that our model atom 
accounts for the Rydberg state deexcitation ladder in an accurate 
way. For comparison reasons, we display also a solar disk-center 
intensity profik^J The observed solar line feature is barely vis- 
ible, which further illustrates the stronger emission from K gi- 
ants. 

The difference between departure coefficients in n — 8 — > 7 tran- 
sitions is of similar magnitude as that between 7 and 6, causing 
comparable emission line strengths. The continuum formation is 
shifted about 0.3 dex outwards when compared to the spectral 
region around 12 yum (the H s opacity increases) but the average 
height of formation for the line-center in the 18.83 yum line in 

Pollux is located at log T500 2.5, slightly further in than the 

12 yum lines. This is because the line is merely the third strongest 
%l — > IV transition with Al = -1. We thus predict that the next 
two Rydberg transitions in the chain (8/ — > Ih located at 18.99 
yum and %k — > 7i at 19.03 /mi) should appear even stronger, how- 
ever these were not covered in our observational setup. 

6. Discussion and conclusions 

The model atom extension and the introduction of collisions 
with neutral hydrogen remove the discrepancy between observed 
and modeled emission in our study. Before the model atom was 

3 Observations from the Kitt Peak solar atlas 
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Fig. 6. Observed and modeled line profiles for the 18.83 fim line 
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[3] The feature to the left in the figure is an OH absorption line. 
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Table 1. Summary of observed magnesium emission line prop- 
erties. 







Normalized 




FWHM 






flux maximum 




[km s" 1 ] 




Star, line [yum] 


12.22 


12.32 18.83 


12.22 


12.32 


18.83 


Pollux 


1.28 


1.37 1.20 


6.8 


7.1 


8.2 


Arcturus 


1.24 


1.23 


8.0 


a 


9.7 


Aldebaran 


1.16 


a 


8.0 


a 




Sun 6 


1.09 


1.15 


5.8 


6.1 





" Emission line blended with IFO absorption line. 
* Modeled solar flux profiles (see text). 
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Table 2. Summary of observed emission line properties around 
12 //m (other elements than magnesium) for Pollux and 
Arcturus. 



Element" 


Wavenumber 6 


Normalized 




Wavelength 


flux maximum 




[cm- 1 ] 


[/an] 


Pollux 


Arcturus 


Si I 


810.360 


12.340 


1.07 


1.08 


Si I 


810.591 


12.337 


1.04 


1.04 


All 


810.704 


12.335 


1.07 


1.06 


Si I 


811.709 


12.320 


1.06 


1.05 


Si I 


813.380 


12.294 


1.06 


1.07 


SiF 


814.273 


12.281 


1.04 




Car* 


814.969 


12.270 


1.05 


1.02 


AIT 


815.375 


12.264 


1.03 




SiF 


815.979 


12.255 


1.03 





" Line identifications based on 
Chang & Noves (1983) and Chang 0984) . 

* From the solar observations by Brault & Noyes (1983). 

' Emission line blended with OH absorption line. 

d Line identification in Chang & Noves (1983) 
stated as 'suspicious'. 



extended, we undertook a number of tests to investigate the large 
discrepancies between observations and models. We therefore 
start our discussion by a short summary of these, with the pur- 
pose to simplify future work. 

The lines are sensitive to the photoionization rates as these af- 
fect the recombination to the Rydberg levels. We thus made an 
ad-hoc increase in the photospheric MARCS radiation field un- 
til the modeled Mg I emission lines matched the observed, but 
found that the required additional mean intensities gave rise to 
surface fluxes that by far exceeded observed ones. We included a 
chromospheric temperature-rise in Arcturus and found it to have 
a negligible impact, with the lines forming in atmospheric layers 
below the temperature minimum. We computed MARCS mod- 
els in spherical geo metry and emp loyed the spherical version of 
MULTI, S-MULTI (lHarperlll994l but differences from plane- 
parallel models were small. No attempt to analyze influences 
from atmospheric inhomogeneities has been made in this work. 
A discussi on about how granulatio n affects the solar lines can be 
found in lRutten & Carlssonl dl994l) . 

6.1. The model atom extension 

All Mg I departure coefficients are shifted upwards (increased) 
by the extension of the model atom, but the upward shift is more 
pronounced in the higher energy levels. It is this change in the 
departure coefficient ratio in the line-forming regions that is suf- 
ficient to cause a significant change in results, i.e., higher emis- 
sion peaks for the larger atom. The enhanced collisional cou- 
pling in the uppermost Mg I levels and to the Mg II ground state 
strengthens the cascading process in rather the same manner as 
the model atom with top-levels n — 9 did when compared to 
one reaching only n — 7 (see Sect. I4.2.TT i. Qualitatively, more re- 
combinations enter at the top-levels, channel down through tran- 
sitions that take part of the Rydberg ladder (see Fig.[TJ and set up 
larger departure coefficient differences between adjacent levels. 
The extension thus has a significant impact on the mid-infrared 
emission lines, whereas the overall character of the Mg I statis- 
tical equilibrium remains. 

To include the energy levels with n > 10 seems especially per- 
tinent when applying the model atom on diluted stellar atmo- 



spheres with low surface gravities (as shown by the large dif- 
ferences in the modeled line profiles of the K giants). These are 
more influenced by radiative transitions, and thus the extension 
ensures that the Rydberg level replenishment from the ion state 
is properly accounted for by including top-levels that are fully 
dominated by collisions. 

6.2. Effects from extra collisions 

Higher rates of collisional excitation and ionization affect the 12 
fim lines in the implemented stellar model atmospheres differ- 
ently. Figs. |3] and [4] show how (by the introduction of collisions 
with neutral hydrogen) the Mgl 12 /mi emission is reduced in 
the Sun, increased in Arcturus and Aldebaran and almost un- 
changed in Pollux. Apparently, a homogeneous increase of col- 
lisional rates actually results in stronger emission in the low sur- 
face gravity atmospheres of Arcturus and Aldebaran (which one 
would perhaps not expect since, generally, collisions act to ther- 
malize lines toward LTE). 

We have analyzed this result by computing additional models 
for the Sun and Arcturus, where we included radiatively allowed 
bound-bound collisions with neutral hydrogen for 1) only the 
three transitions li — > 6h, Ih — > 6g, and Ig — > 6/ and 2) only 
transitions with An, Al — -1 and n > 8 and / > 4. The first test 
imposes stronger collisional coupling only in the n — 7 — > 6 
transitions themselves, whereas the second serves to strengthen 
the collisional bound-bound coupling in transitions in the up- 
permost levels that take part of the Rydberg ladder (see Fig.[TJ, 
while maintaining all other rates. 

We quantify the results of this experiment by the normalized flux 
maximum in the 12.22 /mi line. In the first case, both models 
give lower emission, as expected from the direct thermalizing 
effect on the n = 7 — > 6 transitions. The modeled flux maxi- 
mum above the continuum decreased by ~ 30 % and ~ 25 % for 
the Sun and Arcturus respectively, as compared with the model 
excluding collisions with neutral hydrogen. For the second case, 
however, the emission increased with a similar percentage in the 
solar model, whereas in Arcturus, the modeled flux maximum 
doubled its value. Thus, in a star like Arcturus the enhance- 
ment effect from highly excited lines (see Sect. [3j dominates 
the reduction effect from the higher collisional rates in the line 
transitions themselves, so that when introducing collisions with 
neutral hydrogen homogeneously throughout the model atom, 
the outcome is an emission increase (as seen in Fig. [3]). We 
can understand this by noting that, e.g., the ratio between colli- 
sional and radiative deexcitation rates for the Si — * Ih transition 
(supplying the 12.22 fim line) is Cfi/Rjj ~ 15.0 (Sun) and ~ 0.2 
(Arcturus) in typical line-forming layers when collisions with 
neutral hydrogen are excluded. In principle, this means that the 
contribution from hydrogen is needed in the giants to ensure an 
efficient Rydberg ladder. 

6.3. Observations of Rydberg emission lines around 12 fim 

The observed emission-line flux spectra for Pollux and Arcturus 
in the 12 fim region closely resembles the solar limb intensity 
spectrum, whereas the solar disk-center spectrum lacks emis- 
sion features from other elements than magnesium (see Pig. [5j. 
It is evident that strong emission features from the K giants, 
as compared with solar-type dwarfs, appear for more metallic 
Rydberg lines than magnesium, and the observability of differ- 
ent elements so far follows the same pattern as in the Sun. Future 
observations will tell if this observed trend remains for a larger 
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sample. Model tests with K dwarfs indeed predict lower Mgl 
emission for dwarfs than for giants also within the same spectral 
class (supported as well by observations of the magnetically ac- 
tive K dwarf e Eridani, Richter et al., in preparation). 
For cooler K giants, absorption from water vapor starts to influ- 
ence the 12 jum spectrum and we have detected a blend at 12.32 
pm in Arcturus and Aldebaran. Observations from t his spectral 
regio n in the yet cooler M supergiant Betelguese (Rvd e et al] 
2006) reveal no emission lines above the noise level, which is 
consistent with our modeling of the Mg I lin es in this star (using 
the same stellar parameters as in lRvde et alb . 

6.4. Comparison with other studies 

lUitenbroek & Noyesl (1 19961) modeled the Mgl 12.32 /im line in 
Arcturus and concluded that the computed line was, roughly, half 
as strong as the observed. When using the 'C92' model, we find 
an even larger discrepancy for the 12.22 fim line (see Fig. [3]). 
As already discussed, these longstanding discrepancies between 
observations and models for K giants are removed when using 
our new model atom. We note also that Uitenbro ek & Novesl did 
not detect the water vapor absorption line, which is blended with 
the Mg 1 12.32 fim line in our Fig. [5] The extension of the model 
atom changes the results in the solar case as well (as compared 
with C92 and lRvde et al.l 2004), placing the Mg I lines in higher 
emission, but the former results are recovered by the introduc- 
tion of collisions with neutral hydrogen. 

We have settled here with the rather questioned, albeit standard, 
Drawin recipe for collisions with neutral hydrogen while we 
await results from more rigorous quantum mechanical calcula- 
tions. We have shown the influence from these collisions on the 
formation of the mid-infrared Mg I emission lines, however we 
stress that it has not been an aim of this investigation to put de- 
tailed empirical constraints on their efficiency. Such a task would 
require a larger set of lines, including also other wavelength re- 
gions. Nevertheless, we may still compare our adopted scaling 
factor to the Drawin formula for Mgl collisions with neutral 
hydrogen, Sh = 10~ 3 , with other values from the literature. 
Mentioned in Sect. 14.21 was the exponential decrease resultin g 
in S H = 3 x 10~ 10 for the 12 fim transitions dZhao et al.lll998|) . 
a model which was later ab andoned by the sam e group in fa- 
vor of a constant Sh — 0.05 dGehren et al.l l2004). inferred only 
from optical lines. In a recen t non-LTE abundance stu dy of mag- 
nesium in metal -poor stars dMashonkina et al. ll2008). Sh = 0.1 
is used. The value we find based on the mid-infrared lines is 
more than one order of magnitude lower than the values obtained 
from these two optical studies. Our model atom has not been ap- 
plied to optical lines, however such a combined study should 
be given high priority in future work. We thus conclude that the 
mid-infrared emission lines from near-by giant stars may be suit- 
able diagnostics for testing atomic input data in future non-LTE 
analyses. 

Finally, as we are now able to model and explain the observed 
emission lines for both dwarfs and giants, diagnostic applica- 
tions regarding stellar disk-averaged magnetic fields are possi- 
ble. 
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